Polaron states of electrons over the anisotropic surface of liquid helium 

Sviatoslav S. Sokolov^'^, Antonio Carlos A. Ramos^, and Nelson Studart^ 

""^Departamento de Fisica, Universidade Federal de Sao Carlos, 13565-905 Sao Carlos, Sao Paulo, Brazil 
^ B. I. Verkin Institute for Low Temperature Physics and Engineering, National Academy of Sciences of Ukraine, 61164 

Kharkov, Ukraine 



Abstract 

The energetics and transport properties of the polaron in the anisotropic surface over liquid helium are investi- 
gated. The localization radii and the energy of the ground and excited states are calculated using the variational 
method within the hydrodynamic model of the polaron. In particular, we have considered maximal anisotropy which 
corresponds to the system of electrons in quasi-one-dimensional channels over liquid helium. The polaron binding 
energy is found and the temperature for the polaron formation is estimated below 0.1 K. Solving the hydrodynamic 
equations for fluid velocities, the polaron mobilities along and across the channel are determined. The possibility of 
experimental observation of polarons by measuring the frequency of spectroscopic transitions as well the mobility as 
functions of the holding electric field is addressed. 

PACS numbers: 71.38.-Fi; 73.20.Dx; 73.90.-Ff 



I. INTRODUCTION 



An electron together with its self-induced polarization in a medium forms a quasiparticle which has been named 
polaron. Besides its importance as a standard theoretical model of a fermionic particle coupled to a boson scalar field, 
the polaron has been observed in some physical systems. In particular, there has been great interest in the search 
of polaron states for surface electrons levitated over liquid helium whose properties turn out the system as a good 
candidate for the formation of a polaron of reduced dimensionality. Predicted theoretically long time ago, ||l|,||| the 
surface polaron over liquid helium has been the subject of great amount of experimental work during the last decades. 

Theoretical approaches to investigate the surface polaron over helium are twofold. One is based on the description 
of the dimple state (electron plus the deformation of the isotropic helium surface due to the pressing field) through 
the minimization of the total energy functional of the dimple which leads to a system of coupled equations of motion. 
[Q,HJ^J^ The transportproperties are evaluated in terms of classical hydrodynamical equations since the polaron has 
a high effective mass. The other one involves the concept of a Frohlich-like polaron (a single electron coupled to 
ripplons). |^ The conductivity was calculated in terms of a force-force correlation function within the linear response 
theory. |^ Despite the great difference between the methods, the final results for the structure of the ripplonic polaron 
in both approaches show fair qualitative and quantitative agreement. pO[ | 

In a seminal paper, Shikin and Monarkha Q have determined the ground state of the electron in an isotropic 
surface deformation through the solution of the Schrodinger equation for the electron trapped in the dimple and the 
mechanical equilibrium equation using the Fourier-Bessel transform due to the circular symmetry of the equations. 
The profile of the surface deformation was taken in the harmonic approximation (HA) which allows them to obtain the 
localization length of the electron from the Gaussian wave function, i.e. the polaron radius. Later Monarkha using 
a variational method (VM) was able to obtain the localization length from the minimization of the polaron energy 
calculated with a trial Gaussian wave function. In these works, the influence of external fields and the thickness of 
the helium film on the static and dynamical properties of the polaron was investigated. Marques and Studart |Tl| ] 
have solved in a self-consistent way the Schrodinger equation and the mechanical equilibrium equation obtaining both 
the electron wave function and the profile of the dimple. The comparison of the numerical results with those from 
the HA and VM shows that the VM provides a more exact description of the electron wave function than the HA. 
More recently Farias and Peelers ||l^ have also used the VM for determining both ground and excited polaron states 
taking into account the effect of a positive impurity charge localized in the substrate that supports the helium film. 

Recently, there has been a growing interest in the effects of a corrugated helium surface in the properties of surface 
electrons. One motivation is to make use of suspended helium films to increase electron densities which are limited 
in the case of bulk helium by a surface instability and by the impossibility of obtaining high-mobility electrons on a 
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thin film, in which high densities could be achieved, due to surface roughness of the substrate. The other one is to 
confine surface electrons in one and zero-dimensions as have been realized in semiconductor heterostructures. 

Quasi-one-dimensional (QID) electron systems on the surface of liquid helium have been realized by either geometric 
or electrostatic mechanisms which provide the distortion of the helium surface and a confining electric field holds the 
surface electrons along the formed liquid channels. Multi-wire systems have been created using dielectric optical 
gratings [Q, substrate wrap by nylon threads and metallic gate structures. A single wire was also performed 
using a sharply bent polymer film ]l7t and metallic strips on printed circuit board, jist In such a QID system, the 
electron motion is restricted, in addition to the quantum well due to the holding electric in the direction normal to the 
liquid surface, by a lateral confinement and can be modelled in first order of approximation by a harmonic potential. 
In such conditions, the formation of an asymmetric polaron may become also possible. 

In this paper, we address the question of polaron properties by using the hydrodynamic approach in the case of an 
anisotropic potential U{x,y) due to surface corrugation and in particular for the QID electron system described by 
a parabolic confinement in the y direction. The general formalism can be applied for both circularly symmetric and 
asymmetric polaron states. We consider the properties of both ground and excited polaron states and the anisotropic 
transport properties of the polaron are investigated for surface electrons on He'' and He^ . 

The paper is organized as follows. The general formalism and the main relations are described in Sec. II. In Sec. 
Ill we analyze the properties of the ground and excited polaron states. In Sec. IV we investigate the polaron mobility 
when a driving electric field is applied parallel to the liquid surface. In Sec. V we summarize our main results. 



II. THEORETICAL FORMALISM 



Electrons above the free surface (z = 0) of a helium film with thickness d is prevented to penetrate into the liquid 
because of a high potential barrier (~ 1 eV) at the liquid-vapor interface. If an electric field E± is applied in the 
z direction, the electrons are trapped in a quantum well due to image forces coming from the liquid helium and 
substrate and determined by the potential V{z) — eE±z — Aq/z — Ai/(z -I- d) where Aq — e^(e_ffe — l)/4(e_f/e + 1), 
Ai = e^e//e(es — e_f/e)/(l + effe)^(e_H'e + es) with ene and the dielectric constants of helium and substrate respectively. 
If the barrier height is approximated as infinity, the condition ^{x, y, z = 0) = for the electron wave function must 
be fulfilled for a flat surface. The situation changes drastically if we take into account the surface deformation ^(x, y). 
Now the boundary condition for 'i>{x,y,z) has to be imposed at z = £,{x,y). Shikin and Monarkha |Q show that 
the transformation to a new variable z' = z ~ ^{x,y) allows to avoid the perturbation in the boundary condition for 
^(x,?/, z) leading, however, to modifications in the Schrodinger equation. For ^(0,0) <C d and (z) <C d, where (z) is 
the mean electron distance from the surface, these modifications result, in particular, in the dependence of V{z) on 
an effective holding electric field given by E'^ — E± + A/ed^ with A — e^(es — l)/4(es + 1) where we take ehc — 1, 
and the appearance of an additional term eE^£^{x, y). 

The electron motion along the plane (x, y) in the presence of a magnetic field B in the z direction is described by 
the Schrodinger equation 

1 eB eB 

7^[{Px + -rrvf + iPy - -7r^f\^{x^y) + [<iEli{x,y) + U{x,y)]^{x,y) = eip{x,y). (1) 
2m zc zc 

Here px and py are the x and y components of the momentum operator and we have chosen the symmetric gauge of 
the vector potential A — {—By/2,Bx/2,0). Note that due to the explicit (x, ?;)-dependence of the potential term in 
Eq. (^, Px and py are not conserved. The confinement potential U{x,y) will be considered as a general anisotropic 
parabolic well given by 

U{x,y)^^{ax^+y^). (2) 

where a is the anisotropic parameter. For a = 1, we restore the circular symmetry and a = corresponds the case 
of a QID electron considered in Refs. pT| , p2[ . Then Eq. can be rewritten as 

- I^VV ~ ^L,^ + q {ulx' + coy) + eEl^]i; = e^, (3) 

where V is the 2D-Laplacian, uj^ — uj"^ + 4aci;Q, uiy ~ uj'^ + Aujq with ujc ~ eB/nic the cyclotron frequency, and Lz is 
the angular momentum operator along z which is also not conserved for a ^ 1 when the axial symmetry is lost. As 
a consequence, the second term in Eq. (2) does not contribute to energy eigenvalues e if the ip{x,y) is taken real. If 
we assume that ip is real and vanishes at infinity the electron energy can be calculated from 
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e = 



_2m ^ ■ ■ 8 ■ " ^" 

The total energy of the complex "electron + dimple" can be defined as 



dxdy. 



W 



dxdy, 



(4) 



(5) 



where — pg' /a is the capillary constant, g' — g{\ + 3f /pgd'^), f is the van der Waals coupling constant of the 
helium to the substrate, a and p are the surface tension coefficient and the mass density of helium, respectively, and 
g is the gravitational acceleration. As is seen from Eq. (|^) the inequality > £ is always satisfied. Minimizing the 
Eq. (@), we obtain the mechanical equilibrium equation 



^^(V^e - kU) - eEl^p^ 



(6) 



Note that the quantity of eEVtp'^ plays the role of an electron pressure on the liquid surface. 

As is seen from Eqs. (|3|)-(^, the asymmetry of the electron motion in x and y directions for uj^ ^ i-Oy makes 
inappropriate the use of the Fourier-Bessel transform in polar coordinates as in previous works. Here we use the 2D 
Fourier transform for ^{x,y) as 



i(k^x+kyy) . 



a = a^,y)e-^'^'''^^+''-^^dxdy, 



and the similar transform for ip. Here k is the 2D wave vector and S is the surface area. Using Eq. 
obtain the following expression which connects the Fourier transforms of ^(x, y) and 4''^{x, y): 



eE*, 



(7) 

one can easily 
(8) 



Equations (Q) and (^) can be used to eliminate £,{x,y) from Eqs. ^ and (||) and obtain the total energy of the 
polaron in terms only of the electron wave function and its Fourier transform as 



W 



2cr 



dxdy 



E 



fc2 + P 



(9) 



The electron energy e has the same expression as W except by the absence of the factor 2 in the last term of Eq. (|^). 
For the sake of completeness, we can rewrite the resulting the Schrodinger equation after removing from Eq. 

{§, as: 
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CT \ ^ k +k^ J 

Equation is quite general and can be solved self-consistently and the results are used to evaluate the total 
polaron energy W and surface profile ^{x,y). However this procedure is totally numerical and cumbersome even 
in the symmetric case a = 1. We prefer to make reasonable guesses about the structure of 'ip{x,y) and obtain 
analytical results for the energetics and transport properties of the asymmetric polaron. 

We also hereafter limit ourselves to the most studied case of the QID system model [^,^, where the surface 
electrons are confined by a lateral potential well corresponding the completely anisotropic limit, a = 0, i.e. U{y) = 
toWq7/^/2 with the characteristic frequency defined as luq = a/ eE\_ / niR, where R (~ 10"'' — 10~^ cm) is the curvature 
radius of the liquid in the channel. In the absence of a magnetic field the spectrum for the free electron motion is given 
by E{kx, n) = fi^fc^/2m + fujo{n + 1/2) and the electron localization in the y direction estimated by the parameter 
Lq = (fi/mwo)^^^ <C R and for this reason, curvature effects play no significant role because E"^ pushes the electron 
to the bottom of the channel. This parameter somewhat changes in the presence of a magnetic field. Typical 
values of Lq are of the order of 10~^ cm pl| , p2[ and the condition Lq <C i? is satisfied with great accuracy. Obviously 
for the polaron state the localization parameter along the y direction will be significantly modified. 



3 



III. POLARON ENERGETICS 



A. Ground-state 



Based on the general structure of Eq. ( p^ ) we choose the trial function to describe the ground state of the polaron 



as 



■ exp 



(11) 



where and ly are the electron localization lengths in x and y directions, respectively. Substituting Eq. ( pj] ) into 
Eqs. one obtains the surface deformation ^q{x, y), and the total energy Wq of the polaron ground state as 



£.o{x,y) 



eEl 

47r2(T 



dkii 



exp [- {kill + klil) /4] cos(fc,x) cos{kyy) 



y y 



(12) 



and 



m 



dk^ I dk, 



cxp[- [klll + kltl) /2] 



^2 1 J,2 _i_ J,2 



4m V £2 



(13) 



The ground state energy e'^"^ of electron has the same form except by the coefficient of the integral which in this case 
is (ei?^)2/47r2CT. In the hmit + <C 1, the integral in Eq. (p^ ) can be evaluated analytically resulting in 
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2V2 



(14) 



where 7 = expC, and C = 0.5772... is the Euler-Mascheroni constant. The value of e^^^ is near the same obtaining 
the coefficient 2 instead of 4 in the denominator of first term of Eq. (nj). 

The localization lengths £x and £y, which appear in Eqs. ([ll|)-(p4|)7 have been evaluated by a few methods. In 
particular, the HA similar to that used in Ref. W can be applied and is based on the following approximate expression 
for the surface deformation 



Coix, y) ~ e(0, 0) + - C..(0, 0)x2 + e,jO, 0)y^ 



(15) 



The expansion is similar to that of the potential energy near its minimum value in the 2D oscillatory problem and 
reduces the problem of the electron motion in the dimple to the motion in a parabolic confinement potential, by 
defining and £y as localizations parameters of the 2D harmonic oscillator. We prefer however to use the VM which 
allows to obtain and £y from the minimization conditions of the energy Wq, i.e. dW^/dlx = dWo/diy = of the 
minimum. Note that in the VM the key point for consideration is the expression for the polaron energy, given by 
the Ec^. (|r^), whereas the structure of the dimple potential is less important and is taken into account through Eq. 
(^. From the experimental point of view, using the VM seems more convenient because the characteristic value of 
the helium surface depression in the centre of the dimple is about 10~* cm for ~ 3 kV/cm which is impossible 
to detect directly. Hereafter this estimative is based on actual holding fields in experiments on electrons along QID 
channels on bulk helium, where the effects of the film thickness are neglected. ||lj,|l^ On the other hand, as will see, 
the energy gap between the ground and excited polaron states, calculated using the VM, may in principle be accessible 
experimentally. 

In the VM, and £y can be find from the roots of the system of equations 



1 



1 



1 



and 



1 

^y 



1 



L p^y {^x 



+ iy) 



1 



1 



0, 



(16) 



where = 211(7 /m{eE*^)'^ and L\ — 2h/mu)c- Eqs. ( p6| ) have been solved numerically. The results are presented 
in Fig. 1 as a function of the holding field for some values of the magnetic field and for ^He as the liquid substrate. 
The analytical solution of Eqs. (|l^) is possible in some limiting cases. For B = Q and very high holding fields E*_^ ^ 4 
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kV/cm satisfying the condition Lp <C Lq, one obtains — — V^Lp which is the locahzation parameter for the 
symmetric polaron where the effects of the lateral confinement along y are negligible. For holding fields in the range 
1 < E^l < 3 kV/cm, in the opposite limit Lp ^ Lq, one has £x — Lp and £y ~ Lq which correspond to the numerical 
results shown in Fig. 1. Hence the localization length £y is almost the same as that for the electron moving freely 
along the QID channel. Numerical estimates at _B = are £x ^ 10~^ cm and £y ~ 10~^ cm for 1 < i?^ < 3 kV/cm. 
These values are significantly smaller than the curvature radius R. It means that not only the condition £y <^ R, 
which supports the validity of our QID confinement model, is fulfilled but also the condition £x R is satisfied. The 
application of the magnetic field leads to a decrease of localization lengths in comparison with those for B — 0. For 
high B (lOc ^ ^0)7 one obtains from Eqs. ( p^ ) £x — 6y — \/2Lp for Lp <C Lb and £x — £y — Lb for Lp » Lb- We 
point out that £x and £y, calculated using the HA, give different asymptotic values: ~ ~ Lp for Lp ^ Lq, Lb 
and for _B = 0, the results are £x — Lf/^/2 and £y ~ Lq at L_f » Lo- For high magnetic fields and Lp ^ Lb, one 
has tx — — Lb- It is interesting also to note that if the VM is used for the energy of the electron trapped in the 
dimple e'^"^ instead of Wq, the results are the same as those obtained in HA. One can conclude the solutions of Eqs. 
( p^ ) yielding £x = £y, which is the case of a symmetric polaron, appear either in the limit of very high holding field 
(small Lp) 01 for high magnetic fields {uJc ^ ojq, Lp :$> Lb)- 

Using the localization lengths from Eq. (p^), we depicted in Fig. 2 the polaron energy as a function of the holding 
field for some values of the magnetic field. For £x — Lp and £y ~ Lq and B — Eq. (|lj) is rewritten as 



(eE* )2 , 2^/2 hujo 

In 1 1 

A-Ka y/^kcLp AmL'^p 2 



W^o--^^ln-|^ + ^+^. (17) 



Defining the binding energy E\j as the energy of the polaron state minus the electron energy hLo^j^ of the free motion 
in the lowest subband of the lateral potential, one obtains 



Ana ^J^kcL p AmL 



Note that we have considered very low temperatures T <C hujo where only the lowest subband n = is occupied (for 
the classical SE system ^ \/mT /K). The polaron state is preferable energetically at T < otherwise thermal 
motion can liberate electron from the dimple. We estimate Ef, ~ —0.03 K and E^ ~ —0.3 K for holding electric fields 
of 1 V/cm and 3 kV/cm, respectively, in the case of a ^He substrate. One can conclude that the energetic conditions 
for the formation of the polaron in the QID electron system on the liquid helium surface are almost the same as in 
the case of the 2D symmetric polaron. 

For magnetic fields satisfying the conditions lOc S> 1^0 1 Lp <C Lb, the binding energy, defined by extracting the 
cyclotron energy ?iaJc/2 of the lowest Landau level from Wo in Eq.(^^ can be written as 

47rcr .J^KcLb 

which agrees with the asymptotic value of E^ for the symmetric polaron {loq — 0) for very high B- For B — 5.5 T 
{iOc = 10^^ s~^), we estimate Lb = 1.52 x 10~^ cm and Ei, ~ —0.36 K for E*^ = 3 kV/cm, which is very close to that 
for B = 0. 

We have also estimated the profile of the surface from Eq. (p^). For large distances such that k^x and key ^ 1, 
£,{x,y) decreases exponentially as a function of the distance r = -^/x^ + in the same manner as for symmetric 
polaron. [pl| The value of the dimple depth in its centre can be written as 



-8 



cm 



For B = and 4 ~ Lp, £y ~ Lq, we obtain Co(0,0) ~ -{eEl/2Tra)lTi{A/ ^k^Lp) ~ -1.9 x 10 
for i?^ = 3 kV/cm if the liquid substrate is ''He. For high B, where Lp ^ Lb and ojc ^ t^o, one has 
^0(0,0) ~ -{eEl/2TTa)ln{2/ ^kcLB) ^ -2.2 x 10"* cm. 

It is interesting how to recover from our results those for the symmetric polaron where luq ~ 0. In this case x 
and y directions are equivalent and the z component of the angular moment conserves. Due to this one can replace 
the operator L^ of angular moment in Eqs. (^ and (|lO|) by its eigenvalue Iz which leads to the new contribution 
hujJz/'^ to the electron energy e in the case of non-zero magnetic field. We should put = in the ground state 
of the polaron with lowest angular momentum. Moreover this contribution is zero for real electron wave functions 
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in the polaron state even for nonzero angular momentum. If uq = 0, the solution of Eqs. ( |l6| ) is £ 
1/P — 1/ 4L|, + ^1/16L|^ + and the electron wave function for the ground state is written as 

1 



with 



V'o(a;,y) 



which leads to the same results for Wq, ^oii") as previously found in the case of a symmetric polaron. |^J|,|Tl| In 
particular, the polaron energy is 



Wo 



In 



1 



One can easily also obtain that the second-order derivative of ^(r) at the centre of the dimple [^^"^(0, 0)]q = eE'^ /2tt(t£'^ 
which shows the existence of the minimum at r = 0. 



B. Excited states 



We must choose a trial excited-state wave function orthogonal to ground state wave function given by Eq. (|1] 
which as we have seen is the same as wave function for the 2D asymmetric harmonic oscillator in Cartesian coordinates. 
Hence it is natural to propose the wave function of excited states of the 2D harmonic oscillator as the trial functions 
for the excited polaron states: 



■ exp 









H 







(21) 



Evidently the another excited state |0, 1) can be considered with wave function ^mix^y). The wave function tAqi has 
the same form as V'lo by replacing a; by y in Eq. We find however that the state |1, 0) has energy smaller than 

|0, 1). For this reason the state |1, 0) should be considered the first excited polaron state and we are looking for results 
for this state. The results for the state |0, 1) can be easily obtained in a straightforward way. 

Following the same procedure as for the ground state we arrive to following expressions of the polaron energy and 
the profile of the dimple in excited state: 



m 



dk^ I dk, 



(l-fc2j2/2)^xp [- 



/2] 



1,2 

r\j,y. 



U2 



4to 



3 

si 



1 

"y 



and 



Cio(a;,y) 



dkx I dk 



(1 - kl5ll2) exp [- {ejl + kl5l) I A] cos(fc.x) cos(fc,2/) 



h.2 



1,2 



fc2 



(22) 



(23) 



As before the electron energy e'^^"^ has the coefficient —{eE'^)'^/ATT^a in first term in Eq. (p2|). 

The X dependence of ipioix,y) contributes strongly to decrease the electronic pressure on the liquid surface at 
X = 0, at the centre of the dimple. According to Eq. (^, the pressure is proportional to ipioi^^ v)- As a consequence, 
in the hmit k1{5l -\- 5y) ^ 1, the second-order partial derivative [C"a;]io(0, 0) calculated by Eq. ( p3| ) [^"a;]io(0, 0) = 
— [Cyy]io(0,0) ~ —eE*^/['Ka{5x + 5yY]. This result is a direct consequence of the structure of the trial wave function 
"010(2;, y) and means that the function ^io(a;, y) has a maximum at a; = for fixed y whereas a minimum at y = at 
fixed X and one can not conclude definitely about the nature of the dimple potential as a function of two variables x 
and y ai X = y = Q. In such condition the use of the HA to excited states is inappropriate because it is based on the 
expansion given by Eq. ([l5| ) near the minimum of the dimple potential a.t x — y — Q. For this reason the VM is the 
only consistent way to calculate the polaron properties in the excited state. 

As before the localization lengths 5x and 5y can be obtained, in VM, from the conditions of a minimum of energy 
Wia given by Eq. (H). For kl{5l + (5^) < 1, we obtain 
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5x{25x + 'iSy 



^/lkc{6x + Sy) 



4(4 



Sy)' 



1 

4m V S'i. Sy 



(24) 
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Imposing that dWio /dx = dWio /dy= we arrive to the system of equations 



_3_ 

Si 



L%6x{Sx + Sy) 



1 



Sy L-p5y{5x 



■Sy) 



SxjSx + iSy 
^{Sx + Syf 



= 



^0 



(25) 



In the Hmit of very high holding fields where Lp ^ Lo,Lb, we found 6x — 2.12Lf; Sy ~ 1.73Lf- The localization 
parameters obtained from the minimization of e^^"^ are \/2 times smaller. If S = 0, the analytical solution of Eqs. 
( p5| ) can be fomid in the limit Lq <^ Lp: Sx — V^Lp; Sy ~ Lq. Finally for high B {lUc 3> ujq) one obtains Sx — 
Sy ~ Lb- 

We now discuss the energy gap e*^^"-* — e*^*^-* between the ground and excited states of the electron trapped in the 
dimple. The numerical results of the spectroscopic transition frequency ujio~o — e^^"^ — e^^'^ /h are presented in Fig. 
3 as a function of the holding field for zero magnetic field and two substrates '^He and ^He. In the range 1 < i?^ < 3 
kV/cm, the curves can be described by the analytical expression wio-o — {eE*^Y /2TTah, and changes from 9.8 x 10^ 
s~^ to 8 X 10^ s~^ for the polaron over ^He which corresponds to increasing electron energy from 7.5 x 10^^ K to 
6.2 X 10~^ K. This increase is significantly smaller than \Eb\ calculated by Eq. (|l8|), and the electron transition from 



ground to excited states does not destroy the polaron state. It is interesting to note that Wiq — Wq — [e 



(10) 



.(0) 



]/2 



in such conditions is significantly smaller than \Ei)\ and that also Ci^io-o ^ i^Oi which characteristic values vary from 
5.9 X 10^° s^^ to 1.0 X 10^^ s^^ when E*^ increases from 1 to 3 kV/cm. One should emphasize that the spectroscopic 
frequencies Ci^oi-o for the transition from the ground to the excited state |0, 1) are significantly higher than low-i and 
are the same as loq. It means that after such a spectroscopic transition the electron energy increases from 0.4 K to 0.7 
K when E*^ lies in the range of (1 — 3) kV/cm. These energies are significantly higher than |i?b| and this transition 
should destroy the polaron state. For this reason only the state |1,0) can be considered as the excited polaron state. 
The surface deformation for this state is given by 



eio(0,0) 



eEl 

2-KG 



In- 



^kc{Sx + Sy) S 



(26) 



which follows from Eq. (23) in the limit of k^JS'^ + Sy The absolute values of fio(0,0) are smaller than those 



of Coo(0,0). For El =?> kV/cm, Cio(0,0) ~ -1.67 x 10~« cm for 5 = 0. 

The problem of recovering the symmetric case from the previous results in the case of excites states is more 
complicated than for the ground state. If = 0, the excited state can be described by the trial wave function [[zlJlT 



ipi{r, if) 



( 1 



r exp 



(-rV2A2)e*'-'^ U 



±1, 



(27) 



which is eigenfunction of the angular moment corresponding to the excited state of the 2D harmonic oscillator 
in polar coordinates and A is the single localization parameter. However ^iQ(x,y) is not an eigenfunction of Lz and 
hence cannot be reduced to ipi (r, (p) . For this reason the formal dependence of the excited polaron state on two 
localization lengths Sx and Sy must be kept for luq = 0. In particular Sx 7^ Sy even at very high E^ (Lp Lq) where 
the effects of the potential confinement are negligible and the properties of the symmetric excited polaron state should 
be reproduced. For Sx — 2.12Lp, and Sy ~ 1.73Lf {B = 0), the symmetric polaron energy in the excited state of 
symmetric can be written as 



Wi 



{eEl 



Ana 



In- 



0.550 

krLp 



0.837 



Sy., which means that the influence of B on 



Only at high magnetic fields {Lb ^ Lb and Wc 3> wo) one has Sx 
the electron localization is dominant over confinement effects in the limit when the scale of electron localization is 
determined by the magnetic length Lb- 

It is interesting to know about the structure of the dimple at r = in the limit of loq = 0. To make the problem 



more clear, we consider the second-order derivative 



10 



(0) in the point of r = starting from wave function 



given by Eq. ( pl| ) with Sx — 2.12Lp and Sy ~ 1.73Li?. However for the trial wave function ijjiQ{x,y), [C"r]io depends 
explicitly on ip after conversion to polar coordinates. This dependence has no physical sense for the symmetric polaron 
state and results from the inadequacy of V'lo to describe correctly the excited state in the limit of = 0. Since in 
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the symmetric case, the choice of the coordinate system is arbitrary, one can consider the angle-averaged value of 



the second-order derivative. One can easily show that 



10 



(0)) and ( (0)) are almost zero in the limit of 



10 



kl{6x + Sy) ^ 1. Considering higher-order terms, one finds ([frr^j^o — 3ei?^ /2TTaSlSy which suggests a minimum 
in the centre of the isotropic dimple. This estimative agrees with the results of numerical calculations of Rcf. [^ . 
We have observed that at large distances the function ^i(r) decreases exponentially in the same manner as ^o{r). 

We also can extend our formalism to consider the polaron impurity states treated in Ref . ||l^ . The authors calculated 
the energy gap between the ground and excited electron states over a helium film when a localization potential is 
considered due to a positive impurity charge Ze located on the top of the substrate supporting the film with thickness 
d. It is easy to show that the impurity potential gives a correction Ze/cP to E*^ and an additional parabolic term 
mnj^ {x'^ +y'^)/2 with zu'^ = Ze^/ raS , which turns the problem similar to that considered in the present work. Making 
the necessary adjustments of the expressions of VF^'^^ [s'-^''] and taking the same parameters d,kc and Z as in Ref. 
[ p^ we obtain e'-'^^) — e*^°) ~ 0.60 meV which should be compared with the value 0.445 meV obtained by Farias and 
Peeters in a fully numerical calculation. ||l3| 



IV. POLARON TRANSPORT 



We now investigate the transport properties of the QID polaron. When a driving electric field is applied along the 
plane xy, the surface deformation moves together with the trapped electron inducing a field of hydrodynamic velocities 
in the liquid, which is accompanied by energy dissipation, and leads to a finite value of the polaron mobility. | i|,^, p^ , p4| 
This approach is valid only in the strong coupling limit (high E*j_) where the self-trapped state is energetically favored 
in comparison with the weak coupling limit where the electron is simply scattered by ripplons at T < 1 K. p^ ] In 
order to evaluate the polaron mobility, we employ the energy balance equation eE\\VQ = dps/dt, where vq is the liquid 
velocity at infinity, and pE is the energy density dissipated. The function dps/dt is obtained in a straightforward 
way by finding the normal velocity field induced by the polaron from the solution of the Navier-Stokes equation. The 
calculation procedure in the asymmetric case is similar to that used in Refs. for the symmetric polaron and are 
based in the same system of equations and boundary conditions. However the asymmetry of the surface deformation 
in X and y directions should be taken into account and 2D Fourier transforms like Eqs. (Q) should be employed. If 
the driving field is applied along the x (y) direction, the mobility is given by 



2r^S 



. k 



(28) 



where rj is the helium viscosity. The summation in Eq. (|2^) can be performed analytically and the polaron mobility 
in the case of E\\ along the x direction can be written in terms of the complete elliptic integrals K(t) and E(t) as 



-{- 



[ll - £2)6(4 - iy) 



V2e{Eirv 4 [e (^1 - £l/ei) {iy/lx^K (y^l - q/tl 

-}, 



(29) 



{el - 4)6(4 - 4) 



4 [e ^i-ei/ei) K yi-ei/ei)' 

Here Q(t) is the step function. The general Eq. ( p9| ) is quite simplified in the following limiting cases: 

7r3/2a24 



Px 



and 



V2eiEl)^7] 



7r3/2a24 



if 4 » 4 



Px 



V2e(£;*)2r;ln(4/4) 



if 4 < 4- 



(30) 



(31) 



The polaron mobility along the y direction, py, is the same as for p^ but the index "x" must be replaced by "j/" in 
Eqs. (^)-(|3l|) and vice versa. If ix — iy — i, both px and py reproduce the mobility of the symmetric polaron given 
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As is was seen in Sec. II, the limiting case = f-y = 

can be reached in the hmit of extremely high 
holding fields E\ ^ 4 kV/cm which makes this limit experimentally inaccessible . For 1 < E\ < 3 kV/cm, we have 
£x — Lp ^ (y — Lq. In such a condition the longitudinal the polaron mobility (along the x direction of channel) is 



given by Eq. (30). The transversal mobility /it (along the y direction across the channel) is defined by Eq. (|3l|). As 
it follows from Eqs. ^) and (|l]), ^ii/^t ^ \u{tx/ly) if 4 » iy 

The holding field dependences of /i; and /it are depicted in Fig. 5. The magnetic field influences strongly the 
mobility because it decreases the localization parameters as it can be seen in Fig. 5 where polaron mobilities are 
plotted as functions of E*^ for some values of B. Since ^ ^ and 0-4 ~ 2cr3, the mobilities for the polaron over "^He 
are higher than that for "^He as the substrate. 

V. CONCLUDING REMARKS 

In this work, we have discussed the possibility for the formation of the polaron in the context of anisotropic liquid 
helium surface. In particular, we have evaluated the energetics and transport properties of the polaron in a QID 
channel. We have considered the properties of both ground and excited states of the asymmetric polaron. The 
localization lengths and the ground and excited state energies have been determined as functions of the holding 
electric field within the hydrodynamical model and using a variational approach. We have obtained the frequency 
of spectroscopic transitions between the ground and first levels in the asymmetric polaron states. The transport 
properties of the polaron along the corrugated helium surface have been also studied. 

We think that the measurement of the frequency of spectroscopic transition wio-o between ground and excited 
states (see Fig. 4) should be favorable in experimental attempts to observe the polaron state in QID channel over 
liquid helium surface. Indeed, our estimates indicate that cjiq-o differs significantly of cjq of spectroscopic transitions 
between subbands of free electron state in QID channel which allows to separate the experimental signals for these 
two different types of spectroscopy transitions. 

Another possibility is measuring directly the mobility as in Ref . Q . Our results point out that the mobility of the 
QID polaron is 40% smaller than the mobility of the 2D polaron for the same holding field in the range E*^ < 3 
kV/cm. Our value of the QID polaron mobility ~ 10^ cm^/Vs for T = 0.1 K and E'^ = 3 kV/cm, is three orders of 
magnitude smaller than the mobility of 2D electrons over a flat helium surface and that for electrons moving freely 
along the channel , which allow us to easily distinguish them easily. 

While the above estimates show that our proposed mechanism for the polaron formation in QID systems on helium 
should be compatible with accessible experiments, we add other comments to corroborate or refute our predictions. 
Of particular interest is the possibility to use liquid ^He instead ^He as the substrate for the confined QID electron 
system. [|l^ Recall that for '^He, whose surface tension coefficient is more than twice smaller than that of ^He, the 
binding energy \Eb\ can be substantially larger. Also the spectroscopic transition frequencies are twice larger for '^He 
in comparison with ^He. Furthermore, the description of the polaron in terms of a hydrodynamic viscous model may 
be doubtful for pure ''He for temperatures below 1 K. However, we must emphasize that experimental results [Q have 
been successfully interpreted by theoretical calculations from the hydrodynamic approach for T < 1 K, even though 
this interpretation has been seriously questioned. [|| Furthermore the hydrodynamic approach of viscous phenomena 
is well satisfactory in the case of '^He for T down to ~ 0.1 K. Other attractive substrate should be also the low T 
mixture of "^He and ''He. Recent experiments on QlD-channels have provided a higher effective holding field as 
for electrons on helium films which could greatly enlarging the absolute value of the binding energy E},. p6t 

The electron states have been described in the present work in the one-electron approximation. The validity of 
the results in the case of finite electron densities is limited by the condition that the scale of electron localization 
is small in comparison with the mean interelectronic distance a ~ rij"' along the channel axis where ni is linear 
electron density. For B = and E± < 3 kV/cm where — Lp ^ £y — Lq the characteristic values of £x ^ 10"^ 
cm which conditions rii ^ 10^ cm^'. At higher electron densities, correlation effects could strongly influence the 
polaron properties and the applicability of our theory developed for the low-density limit becomes doubtful. For low 
electron densities, the experimental difficulties may appear in providing the measurements of experimental signal of 
the response of electron system with low level. This signal is proportional to the electron density which is especially 



restricted if we have a single QID channel. One possibility should be to use weakly coupled Nch multi-channels [g^ 16 
like in multiple quantum wells in semiconductor heterostructures. If the average distance between the channels 6 3> a 
one can disregard the correlation between electrons in different channels and the total response should be enhanced 
by the factor Nch- This turns more appropriate the conditions for experimental investigation of QID electron systems 
over liquid helium in the low-density limit. 
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One hope that the asymmetric polaron can be detected at low temperatures, probably at T < 0.1 K, and modern 
experimental methods for achieving low and ultralow temperatures would offer the possibility to observe the polaron 
in the QID channel on the liquid helium surface. 
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FIGURE CAPTIONS 

1/2 

Fig. 1. Localization parameters (a) and (.y (b) in units oi Lp = {2nh^(j/m) /{eE*^) for the polaron ground- 
state on the surface of ^He as a function of the holding electric field E*^ for some values of the magnetic field. 

Fig. 2. Total energy of the polaron ground-state as a function of the holding electric field on the surface of ^He. 

Fig. 3. Frequency of spectroscopic transition from the ground to excited electronic states in the asymmetric dimple 
as a function of the holding electric field for ^He and '^He liquid surfaces. 

Fig. 4. Longitudinal (ni) and transversal (/it) polaron mobilities over '''He as a function of the holding electric field 
calculated using by Eq. (E9^ (solid line) and by the approximate expression [Eqs. d^) and (|l])] (dashed line) for 
zero magnetic field. 

Fig. 5. Longitudinal polaron mobility versus the holding electric field for some values of the magnetic field. 
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